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ABSTRACT 

In this study, the hght curves and spectrum of the photospheric thermal 
radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated using 
2D relativistic hydrodynamic simulations of jets from a collapsar. As the jet 
advances, the density around the head of the jet decreases, and its Lorentz factor 
reaches as high as 200 at the photosphere and 400 inside the photosphere. For an 
on-axis observer, the photosphere appears concave shaped due to the low density 
and high beaming factor of the jet. The luminosity varies because of the abrupt 
change in the position of the photosphere due to the internal structure of the 
jet. Comparing our results with GRB090902B, the flux level of the thermal-like 
component is similar to our model, although the peak energy looks a little bit 
higher (but still within a factor of 2). From the comparison, we estimate that the 
bulk Lorentz factor of GRB090902B is F ~ 2.4 x 102(r/10i2cm) ^here r is the 
radius of the photosphere. The spectrum for an on-axis observer is harder than 
that for an off-axis observer. There is a time lag of a few seconds for high energy 
bands in the light curve. This may be the reason for the delayed onset of GeV 
emission seen in GRB080916C. The spectrum below the peak energy is a power 
law and the index is 2.3 ~ 2.6 which is softer than that of single temperature 
plank distribution but still harder than that of typical value of observed one. 
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INTRODUCTION 



Gamma-ray bursts (GRBs), the most energetic explosion in the universe, show a non- 
thermal spectrum, implying that GRBs originate from an optically thin region (e.g. Piran 
(1999) and references therein). On the other hand, some long GRBs are associated with 
a supernova (SN) explosion (e.g., Woosley and Bloom (2006) and references therein). Are 
GRBs already optically thin when they break out the their progenitors ? The answer is 
obviously no, since the jet density should be high and Lorentz factor should be low. Thus, we 
can expect to observe the thermal radiation from the photosphere prior to GRBs. Even in the 
prompt phase, there may be a thermal radiation component coming from the photosphere. 

About 25 years ago, it was suggested GR Bs should show th ermal spectra unless they are 
coming from relativistically- moving objects (IGoodman Ill986l ). Similar to our motivation, 
the ther mal radiation fr o m th e photosphere o f GRBs were analy ti cally discu ssed in some 
paper s jBlinnikov et al.l Jl999[ ). kvde I fcoosh. IPe'er et all fcoOTh . EH fcoOTh . and iPe'er I 
( 120081 )). In the photospheric model iLazzati et al.l (j2009[ ). thermal radiation is scattered 



by relativistic electrons producing non-thermal photons of GRBs (e.g. Rees and Meszaros 
(1999); loka et al. (2007); Toma et al. (2009 and 2010)). The effects of magnetic fields, for 



examp le th e reconnection, on the emission have been discussed in some papers ( iGiannios 
k00& \ and lzhang fc Pe'er I EA ). 



Thermal components have been observe d in some GRBs f r om th e data obtained from 

J2007h . and iRvde fc Pe'eT ^ 



BATS E catalog and/or BeppoSAX catalog (iGhirlanda et al. 



20061: 



( 2009 ) ) . A thermal comp onent has been reported for GRB060218 associated with SN2006aj f lCampana et al 



Liang et al. II2006I ). Recently, Fermi Gamma-ray Space Telescope o bserved GRB090902B, 



showing a clear spectrum consist of two components f lRyde et al.ll2010l ) with one appearing 
to be a thermal component. 

In this study, we investigate the thermal radiation from the photosphere of GRBs 
through numerical simulations. Some studies have been carried out on the jet propaga- 
tion of GRBs (e.g. Mizuta & Aloy (2009) and references therein); however these studies did 
not determine the photosphere's location and the number of thermal photons radiated from 
GRBs except for Lazzati et al. (2009). In this study, we determine the location and shape 
of the photosphere of GRBs as a function of time and estimate the light curve and spectrum 
of the thermal radiation. Here we present thermal spectrum and light curves for different 
viewing angles and light curves for several energy bins corresponding to Fermi and Swift 
data. We also discuss the unique spectrum of GRB090902B and show the delayed onset of 
hard photons possibly related to the one seen in GRB080916C (lAbdo et al.l 120091 1. and we 
discuss viewing angle effects on short GRBs and X-ray flash (X-ray flare). 
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NUMERICAL SETUP 



This section describes the progenitor model and the numerical conditions of the hydro- 
dynamic calculations. Although the structure of the progenitor at the pre-SN stage is not 



Wooslev & Heeer 


f2006h. 


Yoon & Laneer 


models in 


Wooslev & Heger 


( 


2006h. The 



high spin at the pre -SN s tage; this progenitor model is also used by iMorsony et al.l ( 120071 ) . 
Lazzati et al.l ( l2009l ) and iMorsony et al.l (120101 ). At zero age, the progenitor is assumed to 
have a mass of 16 solar masses and low metallicity (0.01 Zq). At the pre-SN stage, the total 
mass is 13.95 solar masses and the progenitor radius (r*) is 4 x 10^° cm. 

The 2D spherical coordinate system (r x 6) is used for hydrodynamic calculation, as- 
suming axisymmetry and equatorial plane symmetry. The computational domain covers the 
region of 10^ cm < r < 3 x 10^^ cm and 0° < ^ < 90°. We remap the density profile from 
the progenitor model into the computational domain (rmin < f ^ ''"*); assuming spherical 
symmetry. The reflective boundary condition is applied at the polar axis and the equatorial 
plane. Since the stellar wind theory is not well known so far and a weak wind is expected 
due to very low metallicity, we assume the gas outside the progenitor is dilute and the power 
law distribution p = 1.8 x 10~^"^(r/r^.)^g cm""^. Since this weak wind does not affect to the 
jet dynamics and the opacity, the model corresponds to the most luminous case. 



The radial grid consists of Nr = 2640 points, uniformly spaced in logr, extending 



from r„ 



10^ cm to r,. 



3 X 10 cm. The smallest radial grid spacing, at 



IS 



^'"min = 10'' cm, while the largest one, at rmax, is Ar„ 
polar grid is spaced between 0° < 6* < 30°, i.e., A9 = 
grids are spaced in the range of 30° < 9 < 90°. 



= 7.6 X 10^ cm. The uniform 120 
0.25°. The 60 uniform logarithmic 



The mechanism of the jet formation near the central black hole is still under debate (but 
see e.g., Nagataki et al. (2007), Nagataki (2009)). A very hot, relativistic and temporary 
constant jet is injected from the inner most grid, i.e., r = rmin with a half opening angle 
10°. The velocity vector is parallel to the radial unit vector. The luminosity of the injected 
jet is 5.5 X 10^° erg/s. Lorentz factor (Fq) is 5 and specific internal energy (eo/c^) is 80, 
where c is the speed of light. This jet has a potential to be accelerated to a Lorentz factor of 
more than 500, applying relativistic Bernoulli's principle ; hV = const, along a stream line 
(/ioLq = 538 in our case), where h is specific enthalpy, if all thermal energy is converted to 
kinetic energy without dissipation. We follow the jet propagation until the head of the jet 
reaches r = 3 x 10^^ cm. 



A special relativistic hydrodynamic code developed by one of the authors (AM 



Mizuta et al 
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( 120041 ) and iMizuta et al.l fl2006[ )) is used for hydrodynamic simulations. The version with 
2nd order accuracy both in time and space is used, including some dissipation to prevent 
numerical oscillation at strong shocks. 



RESULTS 



3.1. Hydrodynamics 

In the early phase of the evolution, the jet drills through the progen itor envelop e , keep - 



ing good co l 
Zhang et al.l ( 



imation, a s shown in previous numerical simu l ations bv lAlov et al.l ( 120001) . 



2003 



and Lazzati et al 



2004 ). lMizuta et al.l ( 120061 ). iMorsonv et al.l (120071 ). iMizuta fc Alovl f l2009[ ). 
(|2009[ ). The forward shock drives the envelopes to the high pressure and 



high temperature. A cocoon originates from the jet material because of the reverse shock at 
the head of the jet and surrounds the jet. Since the density of the jet is considerably lower 
than that of the st ellar envelopes, a fast backfiow and some vortices appear in the cocoon 
( iMizuta et al.ll2010l ). Internal shocks in the jet appear due to interaction bet ween the jet and 



the hi gh pressure cocoon and backfiows (Komissarov & Falle (1997,1998) and lMorsony et al. 
( I2OO7I )). The jet reaches the progenitor surface at about tiab = 7s. 



When the shock breaks out the progenitor surface, the components near the head of 
the jet starts expanding, resulting in a high Lorentz factor component (F ~ 200) (two top 
panels of Fig. [T| density and Lorentz factor contours at tjab = 30 s). The shocked envelopes 
are also expanding into circumstellar matter. Since most of the components of the jet are 
still surrounded by the expanding cocoon and high density progenitor envelopes, the jet 
remains collimated. The pressure in the cocoon and envelope decreases as they expand. The 
components injected from computational boundary (r = rmin) at a later phase can easily 
propagate in a radial direction without dissipation (four bottom panels of Fig. [T]). The blue 
regions in density contours and the red and purple regions in Lorentz factor contours are 
free expanding regions. The Lorentz factor exceeds 400. This free expansion ends at the 
internal shock. On the other hand, since the velocity of the expanding cocoon and stellar 
envelopes is less than its sound speed (~ 0.5c), it is delayed with the head of the jet. As 
a result, the head of the jet is quite relaxed, resulting in a high Lorentz factor component 
(F ~ 200). The internal structures imprinted in the jet before the shock break still remain 
near the head of the jet. Such discontinuities can be also seen in a ID radial plot (Fig. 2. 
Lazzati et al.l (120091 )). but the structure is different from our simulations. Such difference 
may be caused by the different code and resolution for the hydrodynamic calculations. 
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3.2. Thermal Radiation from Photosphere 

The thermal radiation from the photosphere is derived in post processing, from data 
taken every 1.0 s of the laboratory frame. In each snapshot of the hydrodynamic simulations, 
we find photosphere at unity optical depth for the Thomson scattering. The optical depth 
(r) is defined as follows; 



oo 



crTr(l — /3 COS 9) n dl, 



ph 



where ctt is the Thomson scattering cross section, r(= (1 — is the Lorentz factor, 

/3 is absolute value of velocity normalized by the speed of light, 9 is the angle between the 
velocity vector and the line-of-sight (LOS), and n is proper number density of the electron, 
n = 2p/mHc, where rriHe is the mass of helium atom. Though the progenitor includes many 
elements except hydrogen, we assume that all materials consist of fully ionized helium for 
simplicity. The expression in the integral includes the in verse of beaming factor (6), i.e.. 



r(l — /3cos6')(= 6~^{9)), to include the relativistic effect (lAbramowicz et al.lll99ll ). Here, 
we study the effect of the viewing angles, for an on-axis observer = and for off-axis 
observers, i.e., the angle between the jet axis and LOS, 9^ = 1°, and 2°. 

Assuming an observer located at infinity, the isotropic luminosity of thermal radiation 
from the photosphere is evaluated as 



Liso{9.) = acj 5 {9YT^ cos 9p^dS, 



(2) 



where a is the radiation constant, is t he angle between the LOS and the normal vector 



of the emission surface (iPe'er et al.l 120071 ). The isotropic luminosity in the left hand side 
is considered by the observer time (tobs) which is related with the laboratory time (tiab) as 
^obs = ^lab + d/c, where d is the distance between each photosphere and the observer. The 
local temperature at the photosphere is evaluated as T = (3p/a)^/^, where p is the thermal 
pressure of the fiuid. The photosphere for 9^ = 0° on the plane of the jet axis and the 
observer is shown by the solid lines in Fig. [1] Though we also plot the cases for 9^ = 5°, 
and 10°, we do not show light curves and spectrum for these cases due to too short duration 
at the observer frame (see, sec 3.3). The photospheres for each observer almost overlap at 
^lab = 30 s. As the jet advances, its density decreases. As a result, the photosphere shifts 
further inside the jet (bottom 4 panels in Fig. [1] at tuh = 65 and 90 s). To an on-axis 
observer, the photosphere, in particular appears highly concave. 
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3.3. Light curves 

Hereafter we assume that the burst occurs at the redshift oi z = 1. Figure [2] shows the 
hght curves for the observer at different viewing angles as a function of the observer time. 
The observer detects the radiation at different observer times, even if the radiation comes 
from the same laboratory time, due to the curved photosphere. We set tobs = as the burst 
trigger for each viewing angle. Though we integrate thermal radiation only of the whole 
computational laboratory time (0 < tiab < 100 s) due to limited CPU time, the radiation 
will continue for some cases, especially on-axis case. Thus, we should stress that the light 
curves shown here are not completed yet, though the early phase of the light curve is valid. 
We need to follow longer timescales at laboratory time in the future. 

The light curves exhibit time variability, i.e., for the observers ^ = 0° and 1° second 
(10 s < tob < 22 s) and third (tob > 22 s) phase after the first phases in the light curve 
which exhibits constant l uminosity and conti nues for more than 32 s. The feature of light 



curves is consistent with iLazzati et al.l (120091 ) in which the case of 6' = 0° is shown. Since 



the observer sees the skin of the jet at the beginning, the photosphere moves almost at the 
speed of light. The duration for the observer is about 1/7^ times as long as that for the 
laboratory frame, where 7 is the Lorentz factor of the motion for the photosphere (7 ~ 10). 
The arrival time of the radiation concentrates within a few seconds for the observer which 
causes quick rise. As time passes, the density near head of the jet decreases due to the 
expansion of the jet. The beaming factor for the on-axis observer is very large, since the 
velocity vector is almost parallel to the jet axis. So, the distance between the forward shock 
and photosphere increases, i.e., Eq{T] The observer can see the region very deep inside the jet. 
The photosphere for the on-axis observer is concave as shown in bottom four panels of Fig. 
[U Because of discontinuities in the Lorentz factor and density in the jet, the photosphere 
suddenly moves deeper inside the jet. The second phase caused by the detection of the 
internal shock (F ~ 200) and beaming factor at the photosphere is more than 400. The 
third phase is caused by the radiation from much deeper side. The Lorentz factor at the 
photosphere is smaller than that in second phase. 

Figure [3] shows the photon number flux for an observer in the four energy bands as 
usually used in swift analysis. Figure H] is same as Fig. [3l but different energy bands as 
usually used in fermi analysis. For 6^ = 0°, only soft photons {E < 100 keV) come in the 
beginning (tob < a couple of seconds) since the Lorentz factor of the photosphere is still 
small. Then, high-energy photons follow. A few seconds time lag for high energy photons 
can be clearly seen in Figj4]for both on-axis observer {6^ = 0°) and off-axis observer {6^ = 1°, 
and 2°). 
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3.4. Spectrum 

Figure |5] shows the vF^, plot for an observer {9^ = 0°, 1°, and 2°). The spectrum is 
superposition of Plankian at local rest frame and each of them is boosted by the beaming 
factor. The interval of the whole observer time is integrated (top panel). The distribution 
is power law both below and above the peak energy. For 6^ = 0°, the spectrum has the 
peak energy at about 270 keV. A comparison with Fig. |4] shows that the spectrum is soft in 
the beginning, but becomes hard later. Further, as can be expected from Fig. [3] and HI the 
spectrum becomes softer for larger viewing angle cases {6^ = 1° and 2° cases). 

The spectrum below the peak energy has power law index with 2.3 ~ 2.6 for all viewing 
angles. The indices are softer than that of single temperature plank distribution, see two 
bottom panels in Fig. |5] in which single temperature plank distribution, fitting to the peak 
energy and absolute value, is shown by thin dashed lines. Since the spectrum is superposition 
of different temperature Plank distribution, the power law index is ~ 2.0 near the peak 
energy and it continues to one tenth of the peak energ y. The index is s till harder than 



that of typical value of observed one, i.e.. Band function (IBand et al.l 119931 ). T he index for 



below the peak energy band in the spectrum of GRB090902B is about unity ( iRyde et al. 



20101). If we apply higher angular resolution for hydrodynamic calculation, the spectrum will 
consist of much more different local temperature (T') and beaming factor (6) contributions, 
resulting different observed temperature {ST'). This would make the spectrum softer. On 
the other hand, the index above the peak energy is softer than that of Band function. This 
suggests the non-thermal process is necessary to reproduce the entire GRB spectrum. We 
will investigate this non-thermal component in the near future. 



4. DISCUSSION and SUMMARY 

We have calculated the light curves and spectrum of the photospheric thermal radiation 
from GRB jets using 2D relativistic hydrodynamic simulations. We found that the thermal 
radiation from the photosphere of GRB jets is "observable" even for cosmologically distant 
GRBs. This component can be observed as a precursor or a thermal component in the 
prompt phase. The simulation needs further study to ascertain whether this component can 
be observed in the afterglow phase. 

As for GRB090902B, the flux level of the thermal-like component reaches as high as 
10^ keV cm~^ s~^ at high peak energy (-Epk) (~ 300 ± 100 keV) even though it is a distant 



GRB with z = 1.822 flRvde et al.ll2010r ). If we put our simulated GRB at ^ = 1.822, the flux 



level would be also about 10"^ keV cm ^ s ^ but with lower peak energy Ep]^ ~ 190 keV. To 
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explain the discrepancy between the peak energy of GRB090902B and our model by S{6)T 
(Eq.(2)), this factor should be greater by a factor of 1.57. On the other hand, since the 
flux level is comparable between GRB090902B and our model, it is suggested that the value 
for 5{eYT^dS is comparable (Eq.(2)). Thus we can deduce that {T'^T'^r'^)^^^/ {T'^T'^r'^)MoAei 
will be of the order of unity. Here F is the bulk Lorentz factor of the photosphere and r is the 
radius of the photosphere. This is because 5{6YT^dS can be rewritten as 27r(5(^)^T^r^r~^ 
(beaming factor), and 5{9) ~ F for a face-on observer. Since r ~ 2 x 10^^ cm and F ~ 2 x 10^ 
in our model, we can constrain the value of (r/F) for GRB090902B as 4.1 x 10^. This can 
be written as F ~ 2.4 x W{r/W^cm) in GRB090902B. It is noted that the radius of 
the photosphere depends on many factors in a complicated way such as the power, initial 
Lorentz factor, and mass loading of the jet. In this analysis, the radius of the photosphere 
of GRB090902B is assumed to be of the order of 10^^ cm like the simulation in this study. 

There are many explanations for many GRBs not showing a clear thermal s pectrum. The 



power -law component probably dominates the thermal component, for example iPe'er fc Ryde 
(120101 ). For all viewing angles the derived spectrum below the peak energy is power law with 
the indices 2.3 ~ 2.6. Maybe this thermal component changes to a non-thermal one due to 
the scatterings with non-thermal, high-energy electrons. GRB jets might b e usually domi 



nated by the magnetic fie ld component, suppressing photospheric emission fjZhang fc Pe'er 



20091 : IZhang fc Yanll201ll ). The strict requirement for collimation may also be the reason. 



The time delay of hard photons shown in Figure [3] and H] is of interest, because it 
may be related to the delayed onset observed in some bursts (e.g. GRB080916C, Abdo et 
al. (2009)). It should be noted that very high-energy e mission may not always correlate 



with lower energy ernission , depending how it is created ( iPe'er et al.ll2006l : iGiannios 
Lazzati fc Begelman Il2010l ). 
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Fig. 1. — The density and Lorentz factor contours around the jet axis (z-axis) at tiab = 30 
(top two), 65 (middle two) and 90 s (bottom two). The photospheres for the observer at 
different viewing angles (yellow : = 0, aqua : 6' = 5°, and orange : 9 = 10°) are shown. 
The observer is at infinity and on 2;-axis {6^ = 0) and at the right-top side {6^ > 0). 
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Fig. 2. — The light curves for the observer at different viewing angles, i.e., 9^ = 0° (solid), 1° 
(dashed), and 2° (dotted). It is assumed that the burst occurs at the redshift of 2; = 1. The 
duration of the radiation is longer for the on-axis observers. The light curve for an on-axis 
observer is bright and has time variability. We should stress the light curves still continues 
and the last part of the light curves for each viewing angle are not completed. 
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Fig. 3. — Photon number flux as a function of observer time for each energy bands (from 
top to bottom, 15keV < E < 25keV, 25keV < E < 50keV, 50keV < E < lOOkeV, lOOkeV < 
E < 350keV). It is assumed that the burst occurs at the redshift of z = 1. The cases of 
three different viewing angles 9^ = 0°, 



1°, and 2° are shown. 




Fig. 4. — Same as Fig. [3l but different energy bands 8keV < E < 250keV (top) and 
250keV < E < 50MeV (bottom). 
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Fig. 5. — uFi, plot for different viewing angles 6^ = 0°, 1°, and 2°. Total time averaged 
spectrum are shown (top) , and time averaged plots for two intervals for the on-axis observer 
are shown (middle and bottom) . It is assumed that the GRB occurs at the redshift of 
z = 1. Thin dashed lines shown in middle and bottom panels are single temperature Plankian 
distributions (peak energy and absolute value are fitted.) The power law indices below the 
peak energy are 2.3 ~ 2.6. 



